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Abstract 

An interlayer phase coherence develops spontaneously in the bilayer quantum Hall system at the filling 



C/3 



cd 



factor v = 1. On the other hand, the spin and pseudospin degrees of freedom are entangled coherently in 
the canted antiferromagnetic phase of the bilayer quantum Hall system at the filling factor v = 2. There 
emerges a complex Nambu-Goldstone mode with a linear dispersion in the zero tunneling-interaction limit 
for both cases. Then its phase field provokes a Josephson supercurrent in each layer, which is dissipationless 



T3 
q . as in a superconductor. We study what kind of phase coherence the Nambu-Goldstone mode develops in 



association with the Josephson supercurrent and its effect on the Hall resistance in the bilayer quantum Hall 
£. ' system at v = 1, 2, by employing the Grassmannian formalism. 

o 

q ; PACS numbers: 73.43.-f, 11.30.Qc ,73.43.Qt, 64.70.Tg 

o 
en 



x 



I. INTRODUCTION 



Quantum Hall (QH) effects are remarkable macroscopic quantum phenomena observed in the 
2-dimensional electron system[|li |2fl. They are so special in condensed matter physics that they 
are deeply connected with the fundamental principles of physics. Moreover, QH system provides 
us with an opportunity to enjoy the interplay between condensed matter physics and particle and 
nuclear physics yfl . 

In particular, the physics of the bilayer quantum Hall (QH) system is enormously rich owing 
to the intralayer and interlay er phase coherence controlled by the interplay between the spin and 
the layer (pseudospin) degrees of freedom[|3l |4|]. The interlayer phase coherence is an especially 
intriguing phenomenon in the bilayer QH system [3], where it is enhanced in the limit A S as - ► 0. 
For instance, at the filling factor v — 1 there arises a unique phase, the spin-ferromagnet and 
pseudospin-ferromagnet phase, which has been well studied both theoretically and experimentally. 
One of the most intriguing phenomena is the Josephson tunneling between the two layers predicted 
in Refs.o, [60, whose first experimental indication was obtained in Ref.[7]. Other examples are 
the anomalous behavior of the Hall resistance reported in counterflow experiments QSi |9J] and in 



drag experiments [10]. They are triggered by the Josephson supercurrent within each layer[12J. 
Quite recently, careful experiments 111 ill were performed to explore the condition for the tunneling 
current to be dissipationless. These phenomena are produced by the pseudospins at v — 1, where 
the Nambu-Goldstone (NG) mode describes a pseudospin wave. 

On the other hand, at v = 2 the bilayer QH system has three phases, the spin-ferromagnet 
and pseudospin-singlet phase (abridged as the spin phase), the spin-singlet and pseudospin fer- 



romagnet phase (abridged as the pseudospin phase) and a canted antiferromagnetic phase[|14 



170 (abridged as the CAF phase), depending on the relative strength between the Zeeman en- 
ergy A z and the interlayer tunneling energy A SAS . The pattern of the symmetry breaking is 
SU(4)-HJ(l)<g)SU(2)<g)SU(2), associated with which there appear four complex NG modesil8|]. 
We have recently analyzed the full details of these NG modes in each phase [|19ll. The CAF phase 
is most interesting, where one of the NG modes becomes gapless and has a linear dispersion 
relation (11911 as the tunneling interaction vanishes (A S as — > 0). It is an urgent and intriguing prob- 
lem what kind of phase coherence this NG mode develops. 

In this paper, we investigate the interlayer phase coherence, the associated NG modes, its 
effective Hamiltonian, the Josephson supercurrent provoked by these NG modes and its effect 



to the Hall resistance in the bilayer QH system at v = 1,2, by employing the Grassmannian 



formalism [|18ll . 

The basic field is the Grassmannian field consisting of complex projective (CP 3 ) fields. We 



introduce n CP fields to analyze the v — n bilayer QH system. The CP field emerges when 
composite bosons undergo Bose-Einstein condensationyD. We first make a perturbative analysis 
of the NG modes and reproduce the same results as obtained in Q19|]. We next analyze the nonper- 
turbative phase coherent phenomena developed by the NG mode having linear dispersion, where 
the phase field $(x) is essentially classical and may become very large, which is necessary to 
analyze the associated Josephson supercurrent. We show that it is the entangled spin-pseudospin 
phase coherence in the CAF phase. The Grassmannian formalism provides us with a clear physical 
picture of the spin-pseudospin phase coherence in the CAF phase, and, furthermore, enables us to 
describe nonperturbative phase-coherent phenomena uniformly in the bilayer QH system. 

We then show that the Josephson supercurrent flows within the layer when there is inhomo- 
geneity in $(x). A related topic has been investigated in [20]. The supercurrent in the CAF phase 



leads to the same formula[|12|l for the anomalous Hall resistivity for the counterflow and drag ge- 
ometries as the one at v = 1. What is remarkable is that the total current flowing in the CAF phase 
is a Josephson supercurrent carrying solely spins in the counterflow geometry. We also remark 
that the supercurrent flows both in the balanced and imbalanced systems at v = 1 but only in 
imbalanced systems at v = 2. 

II. THE SU(4) EFFECTIVE HAMILTONIAN 

Electrons in a plane perform cyclotron motion under perpendicular magnetic field B± and 
create Landau levels. The number of flux quanta passing through the system is N$ = B±S/& D , 
where S is the area of the system and $ D = 2nh/e is the flux quantum. There are iV$ Landau 
sites per one Landau level, each of which is associated with one flux quantum and occupies an 
area S/N§ = 2ir£ 2 B , with the magnetic length £b = ^/h/eB±. 

In the bilayer system an electron has two types of index, the spin index (|, \) and the layer index 
(f, b). They can be incorporated in four types of isospin index, a = f|\fJ,,bt,b.L One Landau site 
may contain four electrons. The filling factor is v = N/N$ with N the total number of electrons. 

We explore the physics of electrons confined to the lowest Landau level (LLL), where the elec- 
tron position is specified solely by the guiding center X = (X, Y), whose X and Y components 



are noncommutative, 

[x,Y\ = -%e B . CD 

The equations of motion follow from this noncommutative relation rather than the kinetic term for 
electrons confined within the LLL. In order to derive the effective Hamiltonian, it is convenient to 
represent the noncommutative relation with the use of the Fock states, 

|n) = 4=f( 6t ) n l >> « = 0,1,2,..., 6|0) = 0, (2) 



where b and W are the ladder operators, 

i= vfe ( *- n 6 ^vfc (x+ ' n <3) 

obeying [6, fet] = 1. Although the Fock states correspond to the Landau sites in the symmetric 
gauge, the resulting effective Hamiltonian is independent of the representation we have chosen. 

We expand the electron field operator by a complete set of one-body wave functions ip n (x) = 
(x\n) in the LLL, 

^ a (x) = ^2c a (n)(f n (x), (4) 

n=l 

where c a (n) is the annihilation operator at the Landau site \n) with a = ft,f-J,,bf,b|. The operators 
c a (m) , cJn) satisfy the standard anticommutation relations, 

{c a (m) , c^(n)} = 5 mn 5 a p, {c a (m) , cp(n)} = {c]jjn) , c^n)} = 0. (5) 



The electron field *l> a (x) has four components, and the bilayer system possesses the underlying 
algebra SU(4), having the subalgebra SU sp i n (2) x SU pp i n (2). We denote the three generators of the 
SU sp i n (2) by r^ pin , and those of SU pp i n (2) by rf pm . There remain nine generators T^ pin r^ pm , whose 
explicit form is given in Appendix A. 

All the physical operators required for the description of the system are constructed as bilinear 
combinations of ip(x) and ip^x). They are 16 density operators 

p{x) = ^(x)^(x), S a (x) = ^(x)rr^(x), 
P a {x) = ^(x)r^(x), R ab (x) = l -^{x)rrrl^{x), (6) 

where S a describes the total spin and 2P Z measures the electron-density difference between the 
two layers. The operator R ab transforms as a spin under SU sp i n (2) and as a pseudospin under 

SU ppm (2). 



The kinetic Hamiltonian is quenched, since the kinetic energy is common to all states in the 
LLL. The Coulomb interaction is decomposed into the SU(4) -invariant and SU(4)-noninvariant 
terms 

H c=lj d 2 xd 2 yV + (x - y)p(x)p(y), (7) 

H c = 2 f d 2 xd 2 yV-(x - y)P z (x)P z (y), (8) 

where 

y±( x ) = J-(-L± _ 1 = |, (9) 

with layer separation d. The tunneling and bias terms are summarized into the pseudo-Zeeman 
term. Combining the Zeeman and pseudo-Zeeman terms we have 

#z P z = - f d 2 x(A z S z + A SAS P X + A bias P z ), (10) 

with the Zeeman gap A z , the tunneling gap A S as, and the bias voltage A b i as = eVbi as . 
The total Hamiltonian is 

H = H+ + H c +H ZpZ . (11) 

Note that the SU(4)-noninvariant terms vanish in the limit d, A z , A SA s, Abi as — > 0. 

We project the density operators © to the LLL by substituting the field operator © into them. 
A typical density operator reads 

Rabip) = e-^ 2 /*R ab {p), (12) 

in momentum space, with 

k ^p) = ^ E^i e " pX i m ) ct H r - pmr r n cM, d3) 

ran 

where c(m) is the 4-component vector made of the operators c a (m). 

What are observed experimentally are the classical densities, which are expectation values 
such as p c \p) = (&\p(p)\&), where \&) represents a generic state in the LLL. The Coulomb 



Hamiltonian governing the classical densities are given byQ2JJ]: 

H* = 7T J d 2 P V+(p)p c \-p)p c \p) + Arc J d 2 pV D (p)Pf(-p)Pf(p) 

- \ I d 2 P V*(p)[S: l (-p)S«(p) + Pt{-p)P c a \p) + R a l b (- P )R a l b (p)} 



vr J d 2 P v x (p)[st(-p)st(p) + pf(- P )pf(p) + Ri z {-p)RtAp)\ 

I j ' d 2 pV x (p)p ci (-p)p ci (p), (14) 



where Vd and Vx are the direct and exchange Coulomb potentials, respectively, 

= 2 , V2^e 2 
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e 2 BP 2 /4 



(15) 



with V X = V£ + V x , V& = V$- V x , and 



V±(p) = ^3 (1 ± e^ d ) e -V 2 /2 ) 

„2p2 poo 

/ (£|p 2 /4)e-^ p2 / 4 ± *-p. \ dke-^ k2 - kd J (f B \p\k). (16) 

4 7re J 
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87re 

Here, Io(x) is the modified Bessel function, and Jq(x) is the Bessel function of the first kind. 

Since the exchange interaction V ± (p) is short ranged, it is a good approximation to make 
the derivative expansion, or, equivalently, the momentum expansion. We may set p cl (p) = po, 

Sa(p) = P$"Sa(p), Pa(p) = P^a(p), and B^ h {p) = P^abip) for tne stud y of N G modes. 
Taking the nontrivial lowest-order terms in the derivative expansion, we obtain the SU(4) effective 
Hamiltonian density 

n eff = j d (J2(^ s -f + ( W 2 + {d k n ab f) + 2J- (^{d k s a f + {d k v z f + (d k n a 

+ P4> [e cap (V z ) 2 - 2e x (^(S a ) 2 + {K az f} - (A Z S Z + A SAS V X + A bias P 2 

where p$ = p /u is the density of states, and 

1 



(17) 
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with 



p° - _ 
^ ~ 4ne£ f 



(18) 



(19) 



This Hamiltonian is valid at v = 1,2 and 3. 

It should be noted that all potential terms vanish in the SU(4)-invariant limit, where perturbative 
excitations are gapless. They are the NG modes associated with spontaneous breaking of SU(4) 
symmetry. They get gapped in the actual system, since SU(4) symmetry is explicitly broken. 
Nevertheless, we call them the NG modes. 



III. BILAYER QUANTUM HALL SYSTEM AT v = 1 

In this section, we first show the ground state structure and the associated NG modes. We then 
show the interlayer phase coherence, the associated Josephson supercurrent, and its effect on the 
Hall resistance, in the limit A SAS — > 0. 

A. Ground state structure 

We introduce the CP 3 field based on the composite boson theory. An electron is converted 
into a composite boson by acquiring a flux quantum in the QH state. The CP 3 field emerges 
when composite bosons undergo Bose-Einstein condensation. The dimensionless SU(4) isospin 
densities are given byH: 

S a (x) = intrrn, 



1 



,t^pp in - 



Va{x) = ^n'T^-n, 



(20) 



where n is the CP 3 field of the form n(x) = (n f ^, n f ^, n b ^, n b ^) . 
The ground state at the imbalanced configuration a is given by 



(nf,n*\nf,n^) = (1,0,0,0), 



(21) 



in the bonding-antibonding representation, which reads 
I VTT^ 

o a/TT^ 
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in the layer representation. The ground-state values of the isospin fields are 



l+O-Q \ 

2 



1— ffp 
2 



(22) 



S a — 9 az ' 



V g = - 

a 2 



I - C^az + ^a* ) , K 9 ab = -J>az 



1 - Oo4e + o-q^ , (23) 



all others being zero, giving a unique phase. The residual symmetry keeping the ground state 
invariant is U(3). Thus, the symmetry-breaking pattern is SU(4)— )-U(3). The target space is the 
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FIG. 1 : The lowest Landau level contains four energy levels corresponding to the two layers and the two 
spin states. They are shown in (a) for Asas > Az and (b) for Asas < Az- The lowest-energy level consists 
of up-spin symmetric states in the balanced configurations, and is filled at v = 1. It is the spin-ferromagnet 
and pseudospin-ferromagnet state. Small fluctuations are NG modes r/ s , rj p , and rj T . 



coset space 



CP 3 = SU(4)/U(3) = U(4)/[U(1) g> U(3)], 



(24) 



which is the complex projective (CP) space. 



B. Effective Hamiltonian for the NG modes ati/ = l 

From the previous subsection, we see that the symmetry-breaking pattern is given by (|24l . and 
therefore three complex NG modes emerge, which are described by the CP 3 fields. 

We analyze the perturbative excitations around the ground state. We parameterize the bonding- 
antibonding state as 



n 



B| 



1 ~ ks| 2 - M 2 - I 7 ?!-! 2 , n Bi, = Vsi nAt = Vp> n Al = r] T , 



requiring the commutation relations 



Vi{x),v]{v) = Po 1 &ijK x - y) 



(25) 



(26) 



in order to satisfy the SU(4) algebraic relation. i] s describes the spin wave, r] p the pseudospin wave, 
and i] r the i?-spin wave connecting the ground state to the highest level in the lowest level (Fig. [Q). 
The layer field reads 



/n f A 
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(27) 



\ n^ y 



Expanding 



{n B \n Bi ,n A \n Ai ) = (l,^,^,^) + 
for small fluctuations around the ground state, we obtain 

n* = y / ^ (l " ^s| 2 + |r/ P | 2 + h\ 2 )) + Vp yjl 

(i-l(h\ 2 + h\ 2 + h\ 2 ))-v P 



°"0 f l / 1 + (To . 
, n^=7] s \ h 77 r 



1 - (To 



1 + CTo b i / 1 - 00 



We then set 



Vi{x) 



<Ji(x) +i$i(x) 



(28) 



1 -0Q 

2 ' 

1 + 00 

2 
(29) 



(30) 



where p Q <Ji{x) is the number density excited from the ground state to the ?th level designated by 
(|29l , and #j(cc) is the conjugate phase field, satisfying the commutation relation 



y [<?&), 0j(v)] = ^ijK x ~ !/)■ 



(31) 



We express the isospin field in terms of the CP 3 field (|29l) , 

25a = Us + \ (0p0r + Vr) >^s + ^ Mr ~ Vr) , 1 ~ 2\ Vs \ 2 - 2\ Vr A 

2V a = U ;E (s,p,r),-^ p --(a s ^ r -^ s r ),p 2 (s,p,r)J , 
2^=(^( S ,p,r),- tfr+ I( Ml -V,),M S ,P,r)), 
VRya = (r yx (s,p,r),a r - -(cr s cr v + tf s tf p ),r 2/z (s,p,r) J , 



27e 2a = ( r^(s,p,r),-$p + -(CT s $ r -$ s r ),r 22 (s,p,r) ) , 



(32) 



with 



P*(s,p,r) = y 1 - a 2 - a o^ P - 2y 1 - ct 2 , (|r]p| 2 + |^ r | 2 ) - y(a s cx r + i? s t? r ) 



p«(s,p,r) =cr + Jl- a$cr p - 2<j (|r] p | 2 + |ry r | 2 ) + 



TxxKp) P; r 
r^(s,p,r 
r 22 (s,p,r 

*V*(s,p,r 
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0-Q 
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5- (a s (T p + tf s tf p ) , 



0o 
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er, 
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O"0$s + y 1 - ^O^r - y (^p^r ~ ^p^r) 

l-cr 2 - a a p - 2^1 - a 2 (|r/ p | 2 + |r? s | 2 ) + y (^r + « 



o"o + V 1 - ctq^p - 2o "o (|^p| 2 + |^s| 2 ) - 



(33) 



Substituting these into (fTTT ), we obtain the effective Hamiltonian 



cPkHes 



CL K /Lppin T" J CI K /Lmix ; 



(34) 
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H 



ppin 



:i -^ + a ° J - (^ P ) 2 + f 



4^(i - ^o) + 
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T^al 



<7„ 



+ ^(^ P ) 2 + ^-^* 2 
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7t r 



J+ + a J: 



J+ - ct J: 



[{d k a x f + (^ x ) 2 ] + £ ( A z + ^A SAS ^==^ ) [a\ + tf 2 ] 



1 vT^^ 



- " 2 " - [(<V 2 ) 2 + (^ 2 ) 2 ] + ^ ( A z + ^A SAS 



2 °~vTT^b 
1. VT^~o 



(35) 



2 J ^yr+^ 



[al + tf 2 ] 



-— A SAS (a l a 2 + ^ 1 i3 2 ), 



(36) 



where we change the variables in (pol as 



l + ^o 



-m 



1 - <7 



-?72, Vr 



1 - o-Q / 1 + a 

—^ Vl- \ n V2, 



(37) 



and A bias and e"- 1 are given by 



A hi 



O-Q 



bias 



A 



pAsAS + CTqE V , 



(In 



e cap ~~ 4(e D e x ), 



(38) 
(39) 
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respectively. The pseudospin mode is decoupled from other modes, and from (1351) we have coher- 
ence lengths of the interlayer phase field # p and the imbalanced field a p 



£* ■ =21 

Srmin 



=ppin 



zr 



aA 3 ^ 



n 



A 



SAS 



Snnin ~~ ^ 



B\ 



n[(l-a 2 )J s + a 2 Jf] 



" VV:n ""•'^(I-^ + Asas/v/T^! 



(40) 



The #p mode is gapless for A SAS = 0, though the a p mode is gapful due to the capacitance term 
e u=l . 

^cap • 

On the other hand, from (1361 ) for A S as = 0, the two modes r)\ and t] 2 are decoupled. There exist 
no gapless modes in the Hamiltonian (1361) provided A z 7^ 0. 

C. Effective Hamiltonian for the NG modes in the limit Asas — > 



We concentrate solely on the gapless mode in the limit A S as — > 0, since we are interested 
in the interlayer coherence in this system. We now analyze the nonperturbative phase-coherent 
phenomena, where the phase field ~d(x) is essentially classical and may become very large. We 
parameterize the CP 3 field as 

/ rF(x) \ 



n^ix) 



f e vd ^' 2 ^l + a{x) \ 



1 

V2 



n bt (x) 
\ n H (x) ) 

Then the isospin fields are expressed as 



1 1 

S z [x) = -, V z (x) = K zz (x) = ~cr{x) 











(41) 



V x {x) = TZ zx (x) = -a/1 -a 2 (x)cos'd(x), V y (x) = K zy (x) = --a/1 - a 2 (x) sin#(aj), 

(42) 



with all others being zero. From (|42l we obtain the effective Hamiltonian 

H eff = ^(l - o 2 {x)){d k §{x)f + Uj s + - ^ A {d k a{x)f 



. PoCcap / / \ n2 PoAsAS 

H - A — (o-(x)-cto) o — 



1 - a 2 (a;) 



1 — a 2 (x) cos$(:e) + 



C 



l-a 2 



-a(x) 



(43) 
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The canonical commutation relation is given by 

A) 



} i a(x),$(x)]=i5(x-y). (44) 



From (1431) and (1441 . the Heisenberg equations of motion can be calculated as 

1 



2 2 T d 

hd t #=-d k (J?d k a) + ^a 

Po Po 



(<W " z ~ 2 (d k a) 

1 — a z 



v=1 , acosi!) A o"o 



with 



- e v c -\a - a Q ) - A SAS + A SAS , (45) 

v 1 - cH VI - trg 

/j^ff = --d k {jfd k $) + AsAsv/f^^sin^, (46) 

Po 



2 

jf={i-o*)Ji, j: = j s +-^— 2 J d s . (47) 



D. Josephson supercurrents 

We now study the electric Josephson supercurrent carried by the gapless mode $(cc). In gen- 
eral, the total current consists of three types of current, the Josephson in-plane current J} *, the 
Josephson tunneling current J7~ z Jos , which is proportional to A S as> an d the Hall current l 7 i Ha11 . What 



has been argued in (11311 is that in the case of v = 1, there exists an interlay er voltage Vj unc and 
thus no dissipationless J ] z 0& exists, when a Q 7^ 0. On the other hand, the Josephson in-plane cur- 
rent, which is dissipationless does exist, even for a ^ 0. Here, we assume the sample parameter 
(To 7^ and A S as = so that there is no dissipationless tunneling current J\°* between the two 
layers. 

The electron densities are p} ' = — ep (1/2 ± V z ) = — ep (1 ± cr(x)) /2 on each layer. Tak- 
ing the time derivative and using (1461) we find 

d t p{ = -d t p\ = e -^-WH{x). (48) 

The time derivative of the charge is associated with the current via the continuity equation, 
a t pe (b) = diJ- {b) . We thus identify J- (b) = ±^ Jos (a;)+constant, where 

J!°\x) = ^fdtfix). (49) 

Consequently, the current J} 0ii (x) flows when there exists inhomogeneity in the phase $(x). Such 
a current is precisely the Josephson supercurrent. Indeed, it is a supercurrent because the coherent 
mode exhibits a linear dispersion relation. 

12 



E. Quantum Hall effects 

Let us inject the current J m into the x direction of the bilayer sample, and assume the system 
to be homogeneous in the y direction (Figfl]). This creates the electric field E y <b) so that the Hall 
current flows into the x-direction. A bilayer system consists of the two layers and the volume 
between them. The Coulomb energy in the volume is minimized[12] by the condition E y = Ey. 
We thus impose El = E y = E y . The current is the sum of the Hall current and the Josephson 
current, 

Jl{x) = J jT P ^E y + Jl°\ J h Ax) = ^E y -Jl°\ (50) 



with R K = 2nh/e the von Klitzing constant. We obtain the standard Hall resistance when J^ os = 
0. That is, the emergence of the Josephson supercurrent is detected if the Hall resistance becomes 
anomalous. 

We apply these formulas to analyze the counterflow and drag experiments since they occur 
without tunneling. In the counterflow experiment, the current J m is injected to the front layer and 
extracted from the back layer at the same edge. Since there is no tunneling we have J b = —Jl = 
—J m . Hence, it follows from (|50l) that E y = 0, or 

pf Z7>b 

R xy = ~m = °> R xy — -^ = °- ( 51 ) 

*Jx Jx 

All the input current is carried by the Josephson supercurrent, Jl™ = J7j n . It generates such an 
inhomogeneous phase field that $(x) = (h/eJ®)J m x. 

On the other hand, in the drag experiment, since interlay er- coherent tunneling is absent, no 
current flows on the back layer, or J* = 0. Hence, it follows from (|50l) that Jm = Jl = 
{u/R K )E y , or 

R f = El - ^ (<V\ 

^■xy — ji — y ■• y-J^) 

A part of the input current is carried by the Josephson supercurrent, Jl° s — |(1 — cr )J m . 
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a) cr >0 
A fttttttt 



b) a <0 

* t t t t 



no spin current 



t t T t 



4- 



no spin current 



fttttttt 



FIG. 2: Schematic illustration of the spin supercurrent flowing along the x axis in the counterflow geometry 
for v = 1 bilayer QH system, (a) In the v = 1 bilayer QH system for gq > 0, all spins are polarized into the 
positive z axis. The interlayer phase difference d{x) is created by feeding a charge current J m to the front 
layer, which also drives the spin current. Electrons flow in each layer as indicated by the dotted horizontal 
arrows. The direction of the spin current flowing in the front layer becomes opposite to the direction of 
that flowing in the back layer, and therefore no spin current flows as a whole, (b) In the v = 1 QH bilayer 
system for gq < 0, similar phenomena occur and therefore no spin current flows as a whole. 



F. Spin Josephson supercurrents 

The spin density in each layer is defined by p s ^ m (x) = Sa^ipa, where s c 
and s a — — \h for a — f ^, b \.. By using the formula 



iftfora = f |,b | 



and (l42l) we have 



I Pn( x ) \ 
Pn( x ) 

Pbt( x ) 
\pbi(x) ) 



111 1 1 \ / 

1-1 1 -1 

11-1-1 

\1 -1 -1 1 / 



A) 

2S z [x) 

2P z (x) 

\2R zz {x)J 



( p n {x) \ 

Pn( x ) 

Pb^(x) 
\Pbi(x) J 



A) 
2 



fl + a(x) \ 



l-cr(x) 

\ ) 



(53) 



(54) 
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Then taking the time derivative for p a , we have 



/ d t p$ m {x) \ 



dtp'fl (x) 



spin, 
spin / 



hpo 

4 



( d t a(x) \ 



-d t a(x) 





(55) 



dtpfy ex) 

The time derivative of the spin is associated with the spin current via the continuity equation (in 
this article we neglect the tunneling current): 



dt P s r(x) 



d x Jl^{ 



x 



(56) 



for each a. We thus identify 



j^ix) = -j b s r(x) 



j 



d x ${x) 



all others = 0. 



(57) 



Therefore from (1571 we see that the total spin current j7" Spm = J2 a ^a Pm ls zer ° 5 an d therefore the 
spin Josephson supercurrent does not flow at v = 1 (Fig. El). 



IV. BILAYER QUANTUM HALL SYSTEM AT v = 2 



The standard Hall resistance is given by R^ 



•c.V 



-Rv 



R K at v = 2. On the other hand, it 



has been found experimentally Q8|— llOU that R f = Rk at v = 2. It seems that the interlay er phase 



coherence together with the supercurrent does not develop at v = 2. Note that the experiments 
[8Hl0|] were performed at the balance point a = 0. As we now demonstrate, the interlay er phase 
coherence develops only at the imbalance point a ^ in the CAF phase. 

In this section, we first show the ground state structure and the NG modes for each phase. We 
then discuss the entangled spin-pseudospin phase coherence, the associated Josephson supercur- 
rent and its effect on the Hall resistance in the CAF phase in the limit A SA s — > 0. 
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A. Ground state structure 



It has been shown ll22ll at v = 2 that the order parameters, which are the classical isospin 



densities for the ground state, are given in terms of two parameters a and f3 as 
Si = £(1 - « 2 ) v/T^, VI = *g*a?y/T=P, VI = ^a 2 /3, 

A yy A A 

(58) 

with all others being zero. The parameters a and j3, satisfying |a| < 1 and \/3\ < 1, are determined 
by the variational equations as 



,2 Asas 4e3, (A 2 , - /3 2 A 2 As; 



A2 _ SAS "-A V U r SAS/ /^Q , \ 



A bias _ 4 (e x + 2a 2 (e D - e x )) + 1 



/SAsas A v 7 ^ 3 ^ 2 "' 

where 



A = ^ Al AS a 2 + A|(l - « 2 )(1 - /? 2 ). (61) 

As a physical variable it is more convenient to use the imbalance parameter defined by 

*o = V° z = ^a 2 /3, (62) 

instead of the bias voltage A b i as - This is possible in the pseudospin and CAF phases. The bilayer 
system is balanced at er = 0, while all electrons are in the front layer at a = 1, and in the back 
layer at a = — 1 . 

There are three phases in the bilayer QH system at v = 2. We discuss them in terms of a and 



First, when a = 0, it follows that S° = 1, V° = K Q ah = 0, since A = A ZA /l - (3 2 . Note that 
f3 disappears from all formulas in ((58). This is the spin phase, which is characterized by the fact 
that the isospin is fully polarized into the spin direction with 

S° z = 1, (63) 

all others being zero. The spins in both layers point to the positive z axis due to the Zeeman effect. 

Second, when a = 1, it follows that S z = and (P°) 2 + (V® ) 2 = 1. This is the pseudospin 

phase, which is characterized by the fact that the isospin is fully polarized into the pseudospin 
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direction with 

Vl = VT^, V° Z =P = a , (64) 

all the others being zero. 

For intermediate values of a (0 < a < 1), not only the spin and pseudospin but also some 
components of the residual spin are nonvanishing, and we may control the density imbalance by 
applying a bias voltage as in the pseudospin phase. It follows from (1581 ) that, as the system goes 
away from the spin phase (a = 0), the spins begin to cant coherently and make antiferromagnetic 
correlations between the two layers. Hence it is called the canted antiferromagnetic phase. 

The interlayer phase coherence is an intriguing phenomenon in the bilayer QH system[|3|]. Since 
it is enhanced in the limit A S as — > 0, it is interesting to also investigate the effective Hamiltonian 
in this limit at v = 2. We need to know how the parameters a and (3 are expressed in terms of the 
physical variables. The solutions for (|6TT ) are 



/3 = ±yi-(^f) 2 + 0(A^ AS ), (65) 

with 

A ^ A SA s + 0(A 3 AS ), (66) 

as we shall derive in (|1571 ). By using (1621) we have 

V z =a = ±a 2 + O(Al AS ). (67) 

The parameters a and (3 are simple functions of the physical variables A S as/Az and cr in the limit 
A SAS ->• 0. 

In particular, one of the layers becomes empty in the pseudospin phase and also near the pseu- 
dospin phase boundary in the CAF phase, since we have cr —$■ ±1 as a —$■ 1. On the other hand, 
the bilayer system becomes balanced in the spin phase and also near the spin phase boundary in 
the CAF phase, since we have a — > as a — > 0. 

B. Grassmannian approach 



We employ the Grassmannian formalism [|18ll to make the physical picture of this NG mode 



clearer and to construct a theory which is valid nonperturbatively. The Grassmannian field Z(x) 
consists of two CP 3 fields ri\(x) and n 2 (x) at v = 2, since there are two electrons per one Landau 
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site. Due to the Pauli exclusion principle they should be orthogonal one to another. Hence, we 
require 

nl(x)-n j (x) = 6 ij , (68) 

with i = 1,2. Using a set of two CP 3 fields subject to this normalization condition we introduce a 
4x2 matrix field, the Grassmannian field given by 

Z(x) = (n u n 2 ), (69) 

obeying 

Z ] Z = 1. (70) 

Though we have introduced two fields ^(x) and n 2 (x), we cannot distinguish them quantum 
mechanically since they describe two electrons in the same Landau site. Namely, two fields Z{x) 
and Z'{x) are indistinguishable physically when they are related by a local U(2) transformation 

U(x), 

Z'(x) = Z(x)U(x). (71) 

By identifying these two fields Z{x) and Z'(x), the 4x2 matrix field Z{x) takes values on the 
Grassmann manifold G 4)2 defined by 

C SU(4) n2) 

4,2 U(l) ® SU(2) ® SU(2) " V ; 

The field Z{x) is no longer a set of two independent CP 3 fields. It is a new object, called the 
Grassmannian field, carrying eight real degrees of freedom. 
The dimensionless SU(4) isospin densities are given by 

1 1 2 

S a {x) = -Tr [Z^rrz] = -J2 n \rrn t , 

8=1 

2 



V a (x) = ^Tr [Z^rz] = l -Y^n\rrn h 

i=l 

U ab {x) = ^Tr [Z^rrrz] = i £ 4^^^, (73) 

where n, consists of the basis rii(x) = (rft, rft, n b ^, n b ^) . The ground state is given by Eq. (|58l ). 
which we express in terms of the two CP 3 fields nf. It is straightforward to show that it is given 
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by nf = Unf with 



U = exp 
/ 



T„ 



ppin/ 



+ 



7I\ 



exp 



_ spin ppin/i 
9 a 2/ c 



exp 



spin ppin /i 
r> 'y 'x ^ 



cos fiStHO CQS fc6 

sin ^te) sin ^_£ 
sin fi2fi+=2 CO s &=• 



_sin^te)sin^ 
cos »!E) cos &±^ 



cos fetJE) sin ^ 



-sin^^cos^ cos ^^ sin ^ 
_cos^^sin^ -sin^^cos^ 



!* \ 



cos ^H) cos &±« 



siniHtelsin^ 



ggg+![) ciT1 gg^ga gin Cggg±zO cos e s +e a 



\ — cos v H A ' sin 



where a ,0p, and 6^ are given by 



_sin^±^sin^ cos ^^ cos ^ 



(74) 



cos 



cos^, 



Q = vl — a 2 , sin^o = a, cos Op = \J 1 — /3 2 

A zv ^^^ A 5 . „ _ A SAS 



sin^fl = — (3, 



Ar 



Vl 



a z 



sin^A 



and 



(75) 



(76) 



n? = (1,0,0,0)', n* = (0,0,l,0)'. 

We may introduce perturbative excitation modes iji by introducing the two CP 3 fields rii = Un 
with 



/l-I^f- I^p \ 



ni 



V 

where we parameterize as 



Vi 

V3 



ViW 



( 



n 2 



ivlvi 



~vtv2 



\ 



Vi 



i - \\m\ 2 - \\va\ 2 



V2 



(77) 



V2 



with i — 1,2, 3, 4, obeying the equal-time commutation relations between r^ and 7/j, or 



r]i{x,t),r]Ux,t) 



or 



Po 



2z 



-5y5(aj-y), 



[<Ti(a;,t),^(a;,t)] = —6ij5(x - y). 
Po 



(78) 



(79) 



(80) 



They are required so the SU(4) algebraic relation holds for <S , V a , and S a b. For a detailed discus- 
sion, see Appendix lAl 

We calculate the isospin components (|73l) with the use of n, = t/fij, and substitute them into 
the effective Hamiltonian (fT7l) . In this way we obtain the effective Hamiltonian for rji, which is 
shown to be the same as the one for the NG modes derived in Ref.[ 191. 
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(a) Spin Phase 



(b)Pseudospin Phase 



v=4 



v=3 



r| 



v=2 



v=l 



i 


^ I 


^3 


i 


1 


k 


^2 


A Az 

1 






AsAS 



M 



S{ 



Sf 



v=4 
v=3 

v=2 
v=l 



A A 


i 


k 

i ! 


^2 


A Asas 

^4 






Az 



A\ 



At 



Bf 



FIG. 3: The lowest two energy levels are occupied in the ground state at v = 2. Small fluctuations are 
the NG modes r/i, 772, 773, and 774. (a) For the spin phase, 771 and 772 describe the fluctuation from the 
up-spin symmetric state to the down-spin symmetric state and from the up-spin antisymmetric state to the 
down-spin antisymmetric state, respectively. Their energy levels are degenerated with the Zeeman gap Az- 
On the other hand, 773 and 774, which are fluctuations from the up-spin symmetric state to the down-spin 
antisymmetric state and from the up-spin antisymmetric state to the down-spin symmetric state, have an 
energy gap of Az ± Asas. respectively, (b) For the pseudospin phase 7/1 and 7/2 describe the fluctuation 
from the up-spin bonding state to the up-spin antibonding state and from the down-spin bonding state to 
the down-spin antibonding state, respectively. Their energy levels are degenerated with the tunneling gap 
Asas- On the other hand, r/3 and 7/4, which are fluctuations from the up-spin bonding state to the down-spin 
antibonding state and from down-spin bonding state to the up-spin antibonding state, have an energy gap of 
Asas ± A z , respectively. 

C. NG modes in the spin phase 

As an illustration we study the spin phase at ct = 0, where the transformation (|74T ) is given by 

/lO-l \ 



\/2 



10-1 

10 10 

\0 1 1 / 



(81) 
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by setting a, f3 = 0. We note that 



/ n s ^\ 



n 



n 



S4. 



n 



At 



A| 



/rM 



LTt 



y nr+ J 



n 



fi 



n 



bt 



bi 



iP 



n. 



(82) 



\n°+ ) 



where 



n Sa = -^{n ha + n fa ), n Aa = -^{n ha - n fa ), 



(83) 



with a =|, |. The lowest-energy one-body electron state is the up-spin symmetric state, and the 
second lowest energy state is the up-spin antisymmetric state. They are filled up at v = 2. The 
perturbative excitations r\i are as illustrated in Fig. |3](a). 

It follows from (1731) . (1741) . and (1771) that the isospin densities are explicitly given in terms of 
(Ji{x) and-#j(;r) by 

ai+O-2 _ _ Q _ $1 + $2 _ ~ q ~ _0"i-0- 2 _^ _^l-^2_5 
Ot — t= — = 0"i, O y — j= = 1/1, /<. x:r — ^= = 0"2, /V^x — ;= = #2, 



ft 



:y;y 



V2 " y v/2 

0- 4 - 0-3 _ _ t? 3 - i? 4 _ - 

= — 0-3, /C^ — j= — = 1/3, AC^ 



^2 "^' yx 
0-4 + 0-3 



v/2 



V2 
5. = 1-E 



v/2 



v/2 



= 0-4, 7?. 



^4 + ^3 



! 'M/Z 



v/2 



= 0. 



i- 



1 of + tf? 



i-E 



! 5?+^ 



, ^ = 0-3 0-4 + ^3^4, ^ = 0-4^2 - 0-2^4, 



8=1 i=l 

(84) 

Substituting them into ( TT71 ). we obtain the effective Hamiltonian of the NG modes in terms of the 
canonical sets of er, and #,■ as 



a=l,4 

PoA z 



+ j" E k^) 2 + ( w 



+ 



=1,4 

I O-3O-4 + ^ 3 $4 



+ 



a=2.3 

*£ + /**) E[«2 + £ 



PO^SAS i 



PoA 



bias 



' a=2,3 
52^3 + $2^3 



The annihilation operators are defined by 



tfO*) 



&i{x) +i<di{x) 
V2 



(85) 



(86) 
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with 



„ _ 1/2- 



■&i = P$ 2, &i 



and they satisfy the commutation relations 



[&i(x, t),-dj(y, t)] = iSijSix - y) 



(87) 



(88) 



or 



Vl{x,t),r]f{y,t) =5ij5(x-y), 



(89) 



with i, j = 1,2,3,4. 

The effective Hamiltonian (f8~5l) reads in terms of the creation and annihilation variables (1861) as 



AT a j d 

^° a=l,4 ^° a=2,3 a=l,4 a=2,3 



-St. 



, S T 



+ A bias [r/^ T ^ + r«] - A SAS [^ T r/l + r£ 



»t. 



,«t. 



(90) 



The variables rf 2 , rj^, and ^ are mixing by A S as and A bias . 

In the momentum space, the annihilation and creation operators are tjI k and r]f k together with 
the commutation relations 

\U^l k ]=5 lJ 5{k-k'). (91) 

For the sake of the simplicity we consider the balanced configuration with A bias = in the rest of 
this subsection. Then the Hamiltonian density is given by 



H^ m = / d 2 k -H spin , 



n 



spin 



■2/ spin | ,T/ s P m I -T/ spin 
Tti + 'T-2 ~r ^-3 J 



(92) 



-T/spin _ 


"4 J s ,9 

— fc 2 + A z 
. Po 


sf s 






-2/spin 
rt 2 — 


"4 7 d 

4Js fc 2 + A z + 4e^ 
. Po 


sf s 






-2/ spin 
rt 3 - 


"4 7 d 

s k 2 + A z 
. Po 


+ Ae~ x 


^Jfe^.fe + 


"4 J s l2 
-^fc 2 + A z 

. Po 


st s 
^4,feV4,fe " 


~~ A SAS 



(93) 
(94) 



st s I st s 
V3,fcV4,fe "r ^4,feV3,fe 



(95) 
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We first analyze the dispersion relation and the coherence length of rj\ k . From (1931 ), we have 

(96) 

(97) 



Po 



z- 



e, 



>h 



2/ 



B 



7rJ, 



Ay/ 



The coherent length diverges in the limit A z — > 0. This mode is a pure spin wave since it describes 
the fluctuation of S x and S y as in (|84|) . Indeed, the energy (l96l) . as well as the coherent length (1971) . 
depend only on the Zeeman gap A z and the intralayer stiffness J s . 
We next analyze those of rj\ k : 



EAk) 



S I„2 



u 



'!■> 




k z + A z + 4e 



Xi 



IX 



Ji 



(98) 
(99) 



A z + 4e x 

They depend not only on A z but also on the exchange Coulomb energy e x and the interlayer 
stiffness originating in the interlayer Coulomb interaction. This mode is a i?-spin wave since it 
describes the fluctuation of 1Z XX and 1Z yx . From (l96T ) and (|98T ) we see that, in the one body picture, 
?7^ and r]\ have the same energy gap A z . Indeed, they are described in terms of 771 and 772, having 
the same energy gap A z (Fig. |3](a)). 

We finally analyze those of r^ k and 7/4 k , which are coupled. Hamiltonian (1931 ) can be written, 
in matrix form, 



-2, spin 



Vik 



Vlk 



Ak — A SAS 
-A S as B k 



Vlk 



Vlk 



(100) 



where 



A k 



Po 



A z + 4e 



x- 



Bi 



—k 2 + A z . 

Po 



(101) 



Hamiltonian (1 1001) can be diagonalized as 



^2/ spin 
"•3 




(102) 



where 

E^ = 



Ak + B k 



{A k -B k y + AA 2 



SAS 



^4 



Ak + -Bfe 



(A 



fc - J B fc ) 2 + 4A 2 



SAS 

(103) 
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and the annihilation operator rf ik (i = 3, 4) given by the form 

(V C l + 4A !as + Cfc) ^3,fe - 2A SAS ?Kfc 



^3,fc 



2(^ + 4Al AS + C fcv ^T4A 



2 
SAS 



(V C k + 4A SAS - Cfc) ^3,fc + 2A SA S^4.fe 

'/I/, = V , 7 = , (104) 

; 2 (C 2 k + 4Ai AS - C fcv /^ + 4A^ AS ) 

with Ck = A k — B k . The annihilation operators (11041 ) satisfy the commutation relations 

Vik,vt] =<M(fc-*0, (105) 

with i, j = 3,4. We obtain the dispersions for the modes 77? fc (i = 3, 4) from (1101b and (I103I ). 
By taking the limit fc — > in (11031) . we have two gaps 

Ef =0 = A z + 2e x + [A{e x f + A 2 AS ] * , £? =0 = A z + 2e x - [4(e^) 2 + A 2 AS ] * • (106) 

The gapless condition {E k t^ = 0) implies 

A z (A z + 4e x )-A 2 AS = 0, (107) 

which holds only along the boundary of the spin and CAF phases: see (4.17) in Ref.[|22J]. In the 
interior of the spin phase we have A Z (A Z + 4e x ) — A| AS > 0, which implies that no gapless 
modes arise from 77I and rj\. From (11061 ). in the one body picture, fj\ and fjl have the energy gap 
A z ± A SAS , respectively. Indeed they are described in terms of 773 and 7/4 (Fig. [3] (a)). These 
excitation modes are i?-spin waves coupled with the layer degree of freedom. There emerge four 
complex NG modes, one describing the spin wave (r/l), and the other three the -R-spin waves 

D. NG modes in the pseudospin phase 

For the pseudospin phase, (3 is identified with the imbalanced parameter <To, as we discussed in 
Sect. IIV Al with (l64l) . In this subsection, instead of (3 we express the effective Hamiltonian, the 
dispersions, and the coherence length in terms of cr , since it is a physical variable. 
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From (1741) . by setting a = 1, we have 

/ y/1 + a — s/1 - (T 

1 

U 



1 











\ 











v^^^ v / iT^ 

\ o o -Vi^~o V^T^ J 



(108) 



and 



n 



n 



( n^ \ 

At 



-n 



c/t 



f"4- 



n 



/? 



bt 



t^ 



n. 



(109) 



\n^ ) 



where 



n 



Ba 



1 



(VT^n ba + y/T+^n fa ), 



n 



An 



1 



(vT + ^on bQ - v/T^^on 



to \ 



(110) 



V2 V V U " V2 

with a =|, 4- The lowest-energy one-body electron state is the up-spin bonding state, and the 
second lowest energy state is the down-spin bonding state. They are filled up at v = 2. The 
perturbative excitations r]i are as illustrated in Fig. |3](b). 

We go on to derive the effective Hamiltonian governing these NG modes. From (1731) . (1741 . and 
(1771) . the isospin densities are given in terms of &i(x) and "&%{x) as: 

/ 7 / ,-^ n A -4- ii z \ / 7 / 

^z = OW2 



i--o 2 i-E 



•Pz 



1 - olo 2 + a 1 - Y^ 



a? + tf, 2 



i=l 



i=l 



«->a; 



5"icr 4 + did A ) , <S V = cti^ 3 - ct 3 i?i, 5 2 = 0-30-4 + $ 3 ?9 



■i, 



ft 



■-.'/ 



*i, 



P„=tf 2 



ft™ = 1?H 



-x;y 



ft 



mi 



a 4, 



ft. 



- J I - Org ( CT 2 0- 3 + $ 2 $ 3 + O"0O"3, ^2 = -<T 0- 2 a 3 + # 2 $ 3 ~ a/1 ~ crgflr 3 



ft 



,y.r 



1 - O"0 CT2^4 - 0-4^2 - O"0^4, fty* = CT CT 2 $ 4 ~ Cr 4 2 + \f 1 ~ 0^4 



ft 2 



— i/l - erg cricr 2 + $i$ 2 ) + ctqcti, K zz = -cr a x a 2 + #i#2 ) - V 1 - o"o°"i 



(HI) 



Now, we substitute the isospin densities (11111) into the effective Hamiltonian (TTTl . In this way we 
derive the effective Hamiltonian of the NG modes in terms of the canonical sets of cr, and $j (or 
with (jj and i?j). 
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In the momentum space, this reads 

I d 2 kH p = fd 2 kH p 1 + ! d 2 kH\ + I d 2 ml, 
where 

n\ = A\a{ k cr x>k + B*JlJ hkj 

with Gi fe, and $j >fc given by (f8~71 ), and 



(112) 

(113) 
(114) 
(115) 
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Po 
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SAS 



2a/1 



A4 F 
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;i-*sv. 



°0 J s; 
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-A z /2 
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K -Az/2 
-Az/2 2$ 



\ 



(116) 



We first analyze the dispersions and the coherence lengths from (11141) . since it describes the 
pseudospin wave. It is diagonalized as: 



HI 



d'kElvlUl 



with 



E lk 



Vlk 



Z\JBiCl 



2,fc 



where if 2 k satisfy the commutation relation 



p pt 

V2,fc> 7 / 2 ,fe' 



S(k - k') 



(117) 

(118) 
(119) 

(120) 



Since the ground state is a squeezed coherent state due to the capacitance energy e cap , it is 
more convenient |3|] to use the dispersion and the coherence lengths of &2 and $2 separately. The 
dispersion relations are given by 

AsAS /„ 2 \ nl», 2J 2 A S AS 



7C<) 



E 



O JCTi 

"■2 __ ZZJ_k 2 
Po 



2aA^^ 



+ e cap (l - 0", 



o 2 ), 



^ 



2^,2 

— -fc 2 + 

Po 



2vT 
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0n 
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and their coherence lengths are 



21, 



ttJ? 



\ A SAS _|_ Of , 

\ v^f cap 



c 



21 



B 1 



l *J*y/T=a* 



l-ri) ' "V Asas 

A similar analysis can be adopted for (11131 ), which is diagonalized as: 

Hi = J d*kEXl kV i k 

with 



(122) 



P " ^IBIAI, 



p 

m,k 



i 
7! 



si 



where rj[ k satisfy the commutation relation 



p pt 



Bl 



5(k-k'). 



#i,. 



The dispersion relations of the canonical sets of &i and $i are given by 



K 1 



9 T a ° A 

2J^ k2+ A 



SAS 



Po 2^[\ 

Their coherence lengths are 



2e£(l-*o). 



£ 



CTn 



Po 



*SAS 



2v^^^ 
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(124) 
(125) 

(126) 
(127) 



£ 



*i 



2/ 



B 



ttJJ 7 ' 



^o 2 )' 



T 1 = 2Z 



B1 



^v^ 3 ^ 



7T 



A 



(128) 



SAS 



\ V^l x{ ° 

It appears that ^ CTl is ill-defined for A SAS — > in (1 1281 ). This is not the case due to the relation 
(11301) in the pseudospin phase, which we mention soon. We see that from (11181) and (11241) . in the 
one body picture, rj{ and rfe have the same energy gap A S as- They are described in terms of rji and 
772, having the same energy gap A SAS (Fig. |3](b)). 

Finally, analyzing of the Hamiltonian (11151) as in the case of the spin phase, we obtain the 
condition for the existence of a gapless mode: 



A 



SAS 



T^al 



A 



SAS 



'I -a, 



f" 463,(1 



0"n 



Ai = 0. 



(129) 



This occurs along the pseudospin-canted boundary: see (5.3) and (5.4) in Ref. [22]. Inside the 
pseudospin phase, since we have 



A 



SAS 



v^ 



Un 



A 



SAS 



'1 - ol 



-4e- x (l-*o 



A z>° 5 



(130) 



there are no gapless modes. 
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E. NG modes in the CAF phase 

We derive the effective Hamiltonian of the NG modes in terms of the canonical sets of a* and 
•di. This can be done by substituting (IA8I) and (IA9I ) into the Hamiltonian (TT71) . We first derive the 
Hamiltonian, without taking any limits. Since the expression becomes too extensive, we introduce 
the notation 



c 9a = cos8 a , s 9a = sin8 a , c e = cos 6^, s e = sin Op, c 6s = cos6 5 , s 9s = sm6 s . 



to make the expression for the effective Hamiltonian more manageable. 
Working in the momentum space, the effective Hamiltonian reads 



H c = I d 2 kH c = f d 2 kH\ + f d 2 m 



2 UOJ c 

2> 



where 



AaCn 1 ' 



(131) 



(132) 



2 ,...., ^o'V,;-\,-;- , 12 ,, ,,:.-, M -%S+c^n.i 



u\ = ( ^jffc 2 + A^ *I,A* + -Kj s + si4)k 



cf 



«2 = QHMIQI 



(133) 
(134) 



with 



J? = 4 J s + s^ J s d , M = 4c^e x + A c,-\ 

/ 1? 2 fc \ / A c c c -e c \ 



Q 



hk 
<5"4,fc 



M\ 



C C ^C _JC Q 

-e c -/ c F c 

5 C a c 6 C 

a c £> c d c 

\ b c d c E c J 



(135) 
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The Matrix elements in (11351) are given by 



c c 

and 



2fc 2 

Po 

2fc_ 2 

2tf_ 

Po 
2fc 2 

Po 



2 jp , 2 Td 



M 






2fc 2 
Po 



c e a J3 + S 6» Q J 1 



2c e 



+ 



c G p e a 



P 



~J\ + "^ 2c^e x , 



D c 



2fc' 



Po 



c 6» 5 ( s e Q ^3 + C 6» Q ^l ) + s o s J\ 



Ap %%£<* 



2c e 



13 



c e a J3 + s e a J 1 ) + c e 5 J; 



M 



+ 



M 



,222 r>/ 2 2 , 2 \ 2 

' S 9fs S 9s C e a e cap A \ C 6p S 6 S 1" C flJ S fl_ c 



0« C X' 



^ + 2 , 



2fc « s 2 ^ce s 

c e s c 2e a J 2 "i ^ — e «i ° 



(136) 



r/9 



%^2 + S20B C s e Xi 



Po 

2^ 

Po 

S26 a S29 s l" 2fc2 

4 L Po 



4 



Se s S20 a J2 + L + Cg a S 2 e l3 £on 

Po 4ZA 



JP + JS-JP-J. )+si (2e- x -e. 



-capy 



A' 
~2~ : 



(137) 



with 



^3 — C 8 a J.s + S 9 a Js, J\ — Cg J s + Sg J s , J 2 



— — (J s — J s ), J3 — Ca J s + Sg J s 



S 28 ri 



se s s 2 e a (2e 



X fc cap. 



+ c da - 



A 



SAS 



An 



» e « = 4c L e x + 2s i. e ' 



cap; 



N 



S2e s s 26 a S 6 



~\ X ^cap) 1 



A 



SAS 



A {ce s c da s dp e a + A z ) 



(138) 



2 AA 

where we denote s 2 <j a = sin26 ) Q ,, s 2 ^ = sin 29 p, and s 2 ^ = sin 205. 

It can be verified that the effective Hamiltonian (11331 ) and (11341 ) reproduces the effective Hamil- 
tonian in the spin phase ( 1921) by taking the limit a, ft — > 0. On the other hand, we reproduce the 
effective Hamiltonian in the pseudospin phase (11121) by taking the limit a — > 1 in (|133l) and (|134l) . 

The effective Hamiltonian in the CAF phase is too complicated to make a further analysis. We 
take the limit A S as — > to examine if some simplified formulas are obtained. In particular we 
would like to seek gapless modes. Such gapless modes will play an important role in driving the 
interlayer coherence in the CAF phase. In this limit, we have: 

A SAS 



cos 6^ 



a 



A 7 ' 



°"o 



sin 6^ = ±W 1 



A 



SAS 



A, 



cos 9x = cos 9 n , sin 9x = sin 9 r 



(139) 
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From (1581) and (1139k the classical ground state reads: 



5 



O"0 



, P° = <r , 72°, = sgn(<7 )ft° , ft 



>o 



VKI(i-KI), 



(140) 



all others being zero. We assume a > for definiteness. The transformation (1741 has a simple 
expression: 

(\ 



c/t 







o\ 



(141) 



1 — \CTo\n 



bU 



(142) 



a/1 - |cr | a/K| 

o -vW a/i - kol o 

\0 1/ 

We find n = Wn is of the form (rft, nf \ h *., n^b-h n.^)* by setting 

™fW = (v 71 - Wo\n fi + ^/\a^\n^) 1 nf M = (-^/\a^\n fl 

Consequently, the ground state is such that \rft) and |n^w.) are filled up: The NG modes 771 and 773 
describe an excitation from the state |n f ^) to |rijL^) and \n b ^), respectively, while the NG modes 
7/2 and r/ 4 describe an excitation from the state |n^_ b x) to \n b ^) and |nfi b *}, respectively. A similar 
analysis can be done for a < 0: \n b ^) and |7ifLi ) are filled up, where 

s 



«ftH = (a/ 1 - kol^ ft + \A a o\ n " <■ 



M\ 



S \ f„ LA 



(?o\n K ), 



(143) 



and the gapless mode 774 describes an excitation from the state \n^ h , ) to |njLi ). 

By using ([1391) with (TT331 . and (fT34l) with (fl35l) . (fT36l) . (fTTTl) . and (fT38l we have the Hamil- 
tonian 



4 f 

h = J2 d2kE < 



cf c 



(144) 



together with the dispersion relations (Fig. @]): 

4fc 2 _ . „ 4fc 2 



-El = £?2 



Po 



-Jf + Az, £3 



Po 



J s d + 2A Z + 8 cos 2 9 a e 



x- 



E, 



\k\ 



'8 J? /2fc 2 



(cos 2 2# Q J s d + sin 2 2# a J s ) + 2 sin 2 26 J e 



Po V Po 
where rjf k (i — 1, 2, 3, 4) are the annihilation operators 

$2,fc - «<5"2,fc 



D C X, 



(145) 
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A CT4 



7/lfc 



"J J 2 (0-3,* + M?3,fc) 
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l.fc 
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FIG. 4: Dispersion relations (I145t for the four NG modes J5$. The sample parameters are d = 231, B « 
5.6T, po = 2.7 x 10 15 m~ 2 , and a = 0.1. Inset: Dispersion relations near k = 0. It is clear that E^{k) is 
linear. 



with 



2fc 2 



2ft 2 



A 194 = Jf , A" 4 = (cos 2 26 a J? + sin 2 2fl Q J.) + 2 sin 2 2^(e n - eZ). 



Po Po 

The annihilation operators 77$ *, satisfy the commutation relation, 



'a\ c D C XJ 



c of 



5ij5(k — k') 



1,2,3,4. 



(147) 



(148) 



withi, j 

We summarize the NG modes in the CAF phase in the limit A S as _ > 0. It is to be empha- 
sized that there emerges one gapless mode, rfi k , reflecting the realization of the exact and its 

rp rp 

spontaneous breaking of the U(l) symmetry generated by yx ~ xy . Furthermore, it has the linear 
dispersion relation as in (|1451 ). which leads to a superfluidity associated with this gapless mode. 
All other modes have gaps. 



F. CAF phase in A SA s -^ up to C(A| AS ) 

We focus solely on the gapless mode rj\ ( or 774 ) by neglecting all other gapped modes, and 
derive the effective Hamiltonian for 774 up to 0(A§ AS ). We assume a > for simplicity. 

The two CP 3 fields to be used in the perturbation theory are given by n = Wn with (1741) and 
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W 
Using (l62l) . we can exactly determine as 
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n 2 



1 1 



(149) 
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A 2 
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a 4 + a 2 A 2 (l - a 2 



-a, 



o- 



(150) 



Note that in the limit A SAS — > we obtain (3 — >• 1, which is in accord with our previous calcula- 
tions. Substituting (11501) into ( |59]) . we find 

^2 „,4 i „2A2/1 „,2\ 



A 2 



iSAS n ! + a,Uz< 1 - -'» + ^ " « 4 V^a 2 + A£(l - « 2 ) 



a 2 (a 2 — ctq) 



o J 



/ 



(151) 
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The relation (11511 ) determines the value of a 2 as a function of A z , A SAS , and cr . Substituting 
this value into (11501) we obtain (3 2 as a function of A z , A SAS , cr . We have thus summarized 
our problem into a single equation (I151I ). When A S as is exactly zero, (11511) yields the relation 



a 



| ctq | . Therefore, for weak tunnelings, we search for a solution in the form 



a 



|a | + AA 2 AS + O(A 4 As ; 



(152) 



where we expect A to be a constant. In order to find the value of A we use (11521) and expand the 
relevant combinations in powers of A| AS . In particular, for the first and second terms of (I1511 l we 
find 



AL,« 4 + a 2 A 2 (l-« 2 ) A2 T , (1-AA 2 )A 



^AS 1 - 
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KI)A 2 



a^ 
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8e^A 
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X^Z A 2 

SAS 
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A 2 4Q + 0(A 4 
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Substituting these into (11511) we obtain 



A 2 
^z 
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1 + 



I1-AA 2 )A 2 AS 



A(2-|a |) 



A 2 

^SAS 



- A^^ALs + 0(A 
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;i-|a |)A| |cr |(l — |cr |) |0"o| 

The lowest terms A^ AS disappear automatically. Requiring the Af; AS -terms to vanish, we obtain 

1 Wo 



A 



A z 2(A z + 4e x (l-KD) 



(155) 



32 



and for a 2 we summarize as 

|0b|ll + 



a 



A 2 



2(A z + 4e^(l-H)) A| 



-, ,+0(Al AS ) 
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Using this in (11501) we come to 
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SAS 



A 2 

^z 



+ 0(A 



4 ^ 
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Finally, using (11561) and (11571) in ([75]) and (|60l) . we find: 



sin 2 05 = | cr 1 ( 1 + 



A z + 8e x (l-|cxo|)) A 



SAS 



2(Az + 4e^(l-|a |)) A 2 



0(At AS ) 



Abias = sgn(a )A z 



1 , 4e x + 8(e D -e x )|a | 1 



A 
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A, 



2 A 2 



+ 0(A 



SAS )1 



(158) 
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respectively. Then by using (11561) . (11571 ). (|1581 ). and (11591 ) with ([T71) . we obtain the effective 
Hamiltonian for the gapless mode 7/4 (04 and $4): 



J. 



IAS 



« = %(W 4 ) 2 + ^(Va 4 ) 2 + Ap (e D - e x )\a \ ( 1 - |<x | - ^ S 



J a4 = 2 J s d + 8J s -|a |(l - ko|) + J;(l - 4KI) 
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with 




<7l? 4 = 


/ A 2 

r, 1 jd : t-^SAS 

V Js + Js A 2 



A 



SAS 



A 2 
^z 



(161) 



Taking A| AS = 0, we reproduce the previously calculated expressions (1144 1 and (11451) . 

We wish to derive the effective Hamiltonian for the nonperturbative analysis of the phase field 
•&(x) . For this purpose, it is necessary to start with the parameterization of the Grassmannian field 
valid for arbitrary values of $(x). We make an ansatz 

/ \ / \ 



n 2 
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o ia i9(x) 



J 





3 +t(l-<T O )0(sc) 



(T(X 
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(162) 



/ 



We expand it around -&(x) = and cr(x) = o"o by setting 5a(x) = cr(a3) — cr . Up to the linear 
orders in $(a;) and <5cr(:r), it is straightforward to show that 






(163) 
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where we have set 



(T(X) 



T] 4 (X) 



O"0 



2 y / a (l -a ) 
By requiring the commutation relation (1791 ), we find 

Po 



i$(x)^(l 



O"0 



!j[<r{x),#{v)}=i8{x-v) 



(164) 



(165) 



We have shown that the CP 3 field (11621 ) is reduced to n 2 in (11491 ) in the linear order of the pertur- 
bation fields, apart from the U(l) factor e -' ta ~»v( x ). We may drop it off the parameterization since 
the CP 3 field is defined up to such a U(l) factor. Indeed, such a factor does not contribute to the 
isospin fields. 

Here we parameterize the CP 3 fields as 



ni 



for a(x) > 0, and 



m 



/i\ 

o 
o 

W 

/o\ 
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/ 



n 2 







\ 
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'a(x) 



-M9(as)/2 



y/l-<r{x) 
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(166) 
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n 2 



(e +i W 2 ,/l + a(x)\ 




V e-M^^fMx) J 



(167) 



for o-(jc) < 0. The isospin density fields are expressed in terms of a{x) and -&{x) 

S z (x) = I - \a(x)\, V z (x)=a(x), 

Tl yy (x) = sgn(a )Tl xx (x) = -y/\a(x)\(l - |tr(aj)|) cosi?(ai), 



Ky X (x) = -sgn(a )TZ xy (x) = -y/\a(x)\(l - \a(x)\) sin^(jc) 



(168) 



with all others being zero. The ground-state expectation values are (cr(x)} = a , (i?(x)) = 0, with 
which the order parameters (11401 i are reproduced from (1168) . It is notable that the fluctuations of 
the phase field $(x) affect both the spin and pseudospin components of the i?-spin. This is very 
different from the spin wave in the monolayer QH system or the pseudospin wave in the bilayer 
QH system at v = 1. Hence we call it the entangled spin-pseudospin phase field $(#). 

By substituting (|168l) into (TT71) . apart from irrelevant constant terms, the resulting effective 
Hamiltonianis: 



U 



J,) 



eff 



(W) : 



Jo 



(Vct) 2 + p$e\ ';J-(a-a 



(169) 
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where we have used 



A bias = sgn(a ) [A z + 4e x + 2e^ 1 |cr |] , (170) 

J a = AJ S + p^M- jf, J# = 4J S V |(1 - Icrol). (171) 

|o"o|(l - |o"o|) 

When we require the equal-time commutation relation, 

^[<r(x),#(y)]=i6(x-y), (172) 

the Hamiltonian (11691) is second quantized, and it has the linear dispersion relation 




E, = \k\\l—(—k 2 + 2e» cr A. (173) 

Po \ Po J 

This agrees with E 4 in Eq. (11451) . It should be emphasized that the effective Hamiltonian (11691) is 
valid in all orders of the phase field %9(x). It may be regarded as a classical Hamiltonian as well, 
where (11721 ) should be replaced with the corresponding Poisson bracket. 

The effective Hamiltonian (11691) for i d{x) and a(x) reminds us of the one that governs the 
Josephson effect at v = 1. The main difference is the absence of the tunneling term, which 
implies that there exists no Josephson tunneling. We have shown that the effective Hamiltonian is 
correct up to 0(A| AS ) as A SAS — > 0. Nevertheless, the Josephson supercurrent is present within 
the layer, which is our main issue. 

By using the Hamiltonian (11691) and the commutation relation (11721) . we obtain the equations 
of motion: 

2 T 
Hd t #(x) = — W(aO - 2t:=\a{x) - a ), (174) 

Po 

2 T 

hd t a{x) = — J -±V 2 d{x). (175) 

Po 

G. Josephson supercurrents in the CAF phase 

We now study the electric Josephson supercurrent carried by the gapless mode d(x) in the CAF 
phase, where the further analysis goes in parallel with that given for v = 1. 

The electron densities are p e = — ep (1 ± V z ) /2 = — ep (1 ± cr(x)) /2 on each layer. Tak- 
ing the time derivative and using (11751 ). we find 

dt P { = ~d t p\ = ^V¥(z). (176) 
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The time derivative of the charge is associated with the current via the continuity equation, 

d t p} b) = diJ t m . We thus identify jf b) = ±J7 os (z;)+constant, where 

J^(x) = e -^d^{x). (177) 

Consequently, the current 3^. 0S (x) flows when there exists inhomogeneity in the phase $(cc). Such 
a current is precisely the Josephson supercurrent. It is intriguing that the current does not flow in 
the balanced system since J$ = at a = 0. 

H. Quantum Hall effects in the CAF phase 

Let us inject the current J m into the x direction of the bilayer sample, and assume the system 
to be homogeneous in the y direction (Fig. |5]). By applying the same argument as given in Sect. 
IIII El we show the anomalous Hall resistance behaviours affected by the phase coherence in the 
CAF phase. 

The current for each layer is the sum of the Hall current and the Josephson current, 

J f x {x) = ^E y + Jl°\ J^(x) = ^E y -J^. (178) 

-KK PO -Kk PO 

We apply these formulas to analyze the counterflow and drag experiments without tunneling. With 
the same argument as given in Sect. IIII E[ we have 

E f E b 

R f xy = -£ = o, R b xy = -% = (179) 

J x J X 

in the counterflow experiment. All the input current is carried by the Josephson supercurrent, 
i7i° s = Jm- It generates such an inhomogeneous phase field that 'd(x) = (h/eJ$)J- m x. 
On the other hand, in the drag experiment, we have J m = Jl = (v/R K )E y , or 



El Rk = 1 



R{ = -J- = _E = Rk at v = 2. (180) 



A part of the input current is carried by the Josephson supercurrent, J"J os = |(1 — a )3m- 

In conclusion, we predict the anomalous Hall resistance (11791) and (11801) in the CAF phase at 



v = 2 by carrying out similar experiments [)8l-l 1 0|| due to Kellogg et al. and Tutuc et al. in the 
imbalanced configuration (er ^ 0). 
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FIG. 5: Schematic illustration of the spin supercurrent flowing along the x-axis in the counterflow geome- 
try, (a) All spins are polarized into the positive z axis due to the Zeeman effect at gq = 0. No spin current 
flows, (b) All electrons belong to the front layer at <7o = 1. No spin current flows, (c) In the CAF phase 
for o"o > 0, some up-spin electrons are moved from the back layer to the front layer by flipping spins. An 
NG mode appears associated with this charge-spin transfer. The interlayer phase difference d{x) is created 
by feeding a charge current J- m to the front layer, which also drives the spin current. Electrons flow in 
each layer as indicated by the dotted horizontal arrows, and the spin current flows as indicated by the solid 
horizontal arrow, (d) In the CAF phase for o§ < 0, similar phenomena occur but the direction of the spin 
current becomes opposite. 

I. Spin Josephson supercurrent in the CAF phase 

An intriguing feature of the CAF phase is that the phase field i9(x) describes the entangled 
spin-pseudospin coherence according to the basic formula (11681) . 
Up to 0((cr — o- ) 2 ), we have S z = 1 — |c(a;)|, and we obtain 

«T = dtfjg" = ^[1 + sgn(a )]d 2 J(x), (181) 

dtpt = dtfC = -^l 1 ~ Bgn(a o )]fl£0(aO. (182) 

The time derivative of the spin is associated with the spin current via the continuity equation, 
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d t pf a {x) = d x jf n {x) for each a. We thus identify 



J^(x) = J^(x) = -fd x ti{x), for a > 0, (183) 

J^{x) = j:i u \x) = -y^(x), for <x < 0. (184) 

The spin current J^ m (x) flows along the x axis, when there exists an inhomogeneous phase 
difference §{x). 

In the counterflow experiment, the total charge current along the x axis is zero: Jl{x) + 
J x {x) = 0. Consequently, the input current generates a pure spin current along the x-axis, 

Jr = j** + jjta + j£» + jj«- = sgn(a )^ in . (185) 

This current is dissipationless since the dispersion relation is linear. It is appropriate to call it a 
spin Josephson supercurrent. It is intriguing that the spin current flows in the opposite directions 
for do > and o"o < 0, as illustrated in Figj5] A comment is in order: The spin current only flows 
within the sample, since spins are scattered in the resistor R and spin directions become random 
outside the sample. 

V. CONCLUSION 

In this paper, we have derived the effective Hamiltonian for the NG modes based on the Grass- 
mannian formalism. We have first reproduced the perturbative results on the dispersions and co- 
herence lengths obtained in Ref.[19]. We have then presented the effective theory describing the 
interlayer coherence in the bilayer QH system at u — 1, 2. The Grassmannian formalism shows a 
clear physical picture of the spontaneous development of an interlayer phase coherence. It is to be 
emphasized that the Grassmannian formalism enables us to analyze nonperturbative phase coher- 
ent phenomena such as the Josephson supercurrent. The nonperturbative analysis was beyond the 
scope of Ref.[19]. It has been argued[3] that the interlayer coherence is due to the Bose-Einstein 
condensation of composite bosons, which are single electrons bound to magnetic flux quanta. The 
composite bosons are described by the CP fields, from which the Grassmannian field is composed. 

We have explored the phase-coherent phenomena in the bilayer system. At v = 1, the interlayer 
phase coherence due to the pseudospin, governed by the NG mode describing a pseudospin wave, 
is developed spontaneously. On the other hand, the phase coherence in the CAF phase is the 



38 



entangled spin-pseudospin phase coherence governed by the NG mode $(x) describing the R- 
spin according to the formula (I1681 ). We have predicted the anomalous Hall resistivity in the 
counterflow and drag experiments. It has been shown to exhibit precisely the same behaviour for 
v — 1 and v — 2. The difference between them is that the supercurrent flows both in balanced 
and imbalanced systems at v = 1 but only in imbalanced systems at v = 2. Furthermore, a spin 
Josephson supercurrent flows in the CAF phase in the counterflow geometry, but not for v = 1. In 
other words, the net spin current is nonzero for the CAF phase, while it is zero for v — 1. This is 
due to the spin structure such that the spins are canted coherently and making antiferromagnetic 
correlations between the two layers at v — 2, while the spin is actually frozen and therefore all of 
the spins are pointing to the positive z axis in both layers at v = 1 in the limit A SAS — > 0. 
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Appendix A: Appendix A SU(4) algebra 

The special unitary group SU(N) has (iV 2 — 1) generators. According to the standard notation 
from elementary particle physics[23J, we denote them as X A , A = 1,2,..., N 2 — 1, which are 
represented by Hermitian, traceless, N x N matrices, and normalize them as 

Tr(A A A B ) = 25 AB . (Al) 

They are characterized by 

4 
[X A , X B ] = 2if ABC \ c , {X A , X b } = -TT^d ABC X c , (A2) 

where f A sc and d A Bc are the structure constants of SU(N). We have X A = t a (the Pauli matrix) 
with f A Bc — £abc and d AB c = in the case of SU(2). 

This standard representation is not convenient for our purpose because the spin group is 
SU(2) x SU(2) in the bilayer electron system with the four-component electron field as \I/ = 
(ijF, tp il , ?/> bt , tp bl ). Embedding SU(2) x SU(2) into SU(4) we define the spin matrix by 

^ Pin = ( ^ ° ) , (A3) 

where a = x,y, z, and the pseudospin matrices by, 

( 1 2 \ ( -il 2 \ ( 1 2 \ 

1 t pp- = 1 t pp- = 1 (A4) 

\U o J yn 2 o J \ o -i 2y / 

where 1 2 is the unit matrix in two dimensions. Nine remaining matrices are simple products of the 
spin and pseudospin matrices: 

, T a \ ( ~lT a \ . . ( T a , 

rT n rr= , rrr^ m = , iTifT= ■ (AS) 

Ta \ iT a / ~T a 
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We denote them T a0 = |r^ pin , T 0a = irf to , T ab = |rf in r 6 ppin . They satisfy the normalization 
condition 

Tr(T M!/ T 7(5 ) = 5^5 uS , (A6) 

and the commutation relations 

where f^^s^iv' is the SU(4) structure constant in the basis (IA3I )- (IA5I ). Greek indices run over 
0,x,y,z. 

From (1741) . (1751) . and (1771) . the explicit form of the isospin densities in terms of r\i is given by: 

X 0a; = — cos 9 a sin Oj£Eq x + cos 9 a cos 9/3 cos ^^Xq^ — sin 9 a cos ^ cos 0$T C XX — sin Q sin 9pX c xz 

— cos #„ cos 9/3 sin ^X^ + sin 9 a cos ^ sin 9sll , 
X y = cos 9 s 1q v + sin9 s Z° z , 
X02 = — cos 9 a cos 9pT^ x — cos # Q sin 9p cos ^^X^ + sin 9 a sin ^ cos 9$T XX — sin # a cos 9pX xz 

+ cos # a sin 6p sin 0<s2^ — sin 9 a sin ^ sin 9sXl , 
l x0 = cos9 s T x0 - sin^X^, 
Xrx = — sin # a cos 9^Lq x — sin # Q sin 9$ cos #5Xq 2 — cos 9 a sin 6^ cos 9&T C XX + cos ^ a cos 9pX c xz 

+ sin # a sin 9p sin 0<s2^ + cos 9 a sin 6^ sin ^X^q, 
X = X c 
X X2 = sin 9 a sin fl^gX^ — sin 9 a cos 6^ cos 9§Xq z — cos # Q cos 6p cos ^X^ — cos 9 a sin 0^2^ 

+ sin 9 a cos 6*/3 sin 9sX yy + cos # Q cos 6p sin 9$X c zQ) 
Xyo = cos9 a X y0 - sin 9 a l c zy , 

X yx = - cos 9p sin ^X^ - sin 9pX yx + cos 9p cos ^X^, 
X ra = cos 9 a sin 6> 5 Xq z — sin 9 a sin ^X^ + cos 9 a cos ^X^ — sin 9 a cos ^X^, 
X^ = sin 9p sin ^X^ - cos 9 p X yx - sin 9p cos 6» 5 X^ 2 , 

X 20 = sin 9 a sin ^X^ + cos 9 a sin ^X^ + sin 9 a cos 0a2^ + cos 9 a cos ^X^,, 
X zx = — sin 9p sin 9sX x0 — sin ^ cos 9 5 X C ZX + cos ^X^, 
X 2J/ = sin0 a X£ o + cos 9 a l c zy , 
Z zz = — cos 6p sin 9 s X x0 — cos 6^ cos 9sX c zx — sin ^X^ z . (A8) 
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where we defined l a0 = <S a , X 0a = V a , lab = H ah and 



2S, = Re 



vim + vlv2 - vlvi - viva 



' X 02/ 



Im 



vim + vim - vhi - 772*73 



J xo = Re [ ? 7i+ r 72], X^ = Re [7/3 + 7/4], X^ = Im[77 3 - 7/4], 2£, = Refax - %], 
x lo = Im [*7i + V2} , Z° yx = Im[*7s + ?74] , Zjjv = - Re ^3 - *74] , X£ 2 = Im^i - t/ 2 ] , 



X 1 



:.() 



zy 



i-Ei" 



i=l 



,t. 



* I ' -'-zx 



-Re 



v\m + *?i*72 + Ai + vlv3 



,t„ 



Im 77JT73 + 77^772 + 774*71 + *72*73 



X c 



1*72 1 



1*71 1 



(A9) 



From (|A8I ). (1A9K and the equal-time commutation relations (1791) , it can be verified that the SU(4) 
algebraic relation 



[lfj, u (x,t),Xys(x,t)] = i5(x - y)fv,v,j5,n'v'ZiJ.'v'(y,t), 



(A10) 



is held. 
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